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Abstract. The long term stratospheric impacts due to emis- 
sions of CO 2 , CH 4 , N 2 O, and ozone depleting substances 
(ODSs) are investigated using an updated version of the God- 
dard two-dimensional (2D) model. Perturbation simulations 
with the ODSs, CO 2 , CH 4 , and N 2 O varied individually are 
performed to isolate the relative roles of these gases in driv- 
ing stratospheric changes over the 1850-2100 time period. 
We also show comparisons with observations and the God- « 
dard Earth Observing System chemistry-climate model sim- 
ulations for the time period 1970-2100 to illustrate that the 
2D model captures the basic processes responsible for long- 
term stratospheric change. 

The 2D simulations indicate that prior to 1940, the 
ozone increases due to CO 2 and CH 4 loading outpace the 
ozone losses due to increasing N 2 O and carbon tetrachloride 
(CCI 4 ) emissions, so that ozone reaches a broad maximum 
during the 1920s- 1930s. This proceeds the significant ozone 
depletion during 1 960-2050 driven by the ODS loading, so 
During the latter half of the 2 1 ** century as ODS emissions 
diminish, CO 2 , N 2 O, and CH 4 loading will all have signif- 
icant impacts on global total ozone based on the IPCC AIB 
(medium) scenario, with CO 2 having the largest individual 
effect. Sensitivity tests illustrate that due to the strong chem- 55 
ical interaction between methane and chlorine, the CH 4 im- 
pact on total ozone becomes significantly more positive with 
larger ODS loading. The model simulations also show that 
changes in stratospheric temperature, Brewer-Dobson circu- 
lation (BDC),andage of air during 1850-2100 are controlled «o 
mainly by the CO 2 and ODS loading. The simulated ac- 
celeration of the BDC causes the age of air to decrease by 
~1 year from 1860-2100. The corresponding photochemi- 
cal lifetimes of N 2 O, CFCl;j, CF 2 CI 2 , and CCI 4 decrease 
by 1 1-13% during 1960-2100 due to the acceleration of the es 
BDC, with much smaller lifetime changes (<4%) caused by 
changes in the photochemical loss rates. 
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1 Introduction 

Changes in the atmospheric abundance of halogenated ozone 
depleting substances (ODSs) and the greenhouse gases 
(GHGs) CO 2 , CH 4 , and N 2 O have been shown to sig- 
nificantly impact the chemical and dynamical structure of 
the stratosphere [e.g.. World Meteorological Organization 
(WMO), 2007, 2011]. For example, much of the decline 
of stratospheric ozone during the 1980s and 1990s has been 
attributed to increased atmospheric halogen loading due to 
anthropogenic forcings. Increases in N 2 O and the odd nitro- 
gen species decrease ozone in the middle stratosphere [e.g., 
Crutzen, 1976], while increases in CO 2 and the subsequent 
cooling reduce the temperature dependent ozone loss rates 
and cause ozone increases in the upper stratosphere fe.g., 
Haigh and Pyle, 1979; Rosenfleld et al., 2002]. 

Recent observational studies have detected the beginning 
of the ozone recovery process in the upper stratosphere, 
where ozone is most sensitive to changes in halogen load- 
ing [e.g., Reinsel, 2002; Newchurch et al., 2003]. However, 
detection of the change in the halogen influence on ozone 
can be complicated by the impacts due to long term changes 
in GHGs. It is therefore of interest to separate the relative 
impacts of the different chemical processes that control long 
term ozone changes. 

Another important aspect of the changing atmospheric 
composition impact on the stratosphere is the quantification 
of photochemical lifetimes of the ODSs and GHGs. These 
lifetimes have come under recent scrutiny [Douglass et al., 
2008], as they are important for deriving surface mixing ra- 
tio boundary conditions from emissions estimates for use in 
atmospheric models [Kaye et al., 1994; WMO, 2011]. The 
potential influence on lifetimes of the Brewer-Dobson cir- 
culation (BDC) acceleration due to climate change has also 
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70 been investigated [Butchart and Scaife, 2001 ; Douglass et al., 
2008], 

Most studies of past and future stratospheric change now 
utilize three-dimensional (3D) coupled chemistry-climate 
models (CCMs) [e.g., Eyring et al., 2006; Eyring et al., 

7 s 2007; WMO, 2011]. Some 3D CCM investigations have 
shown the impact of different processes on long term strato- 
spheric change, such as that due to the multi-decadal changes 
in ODS and GHG concentrations and sea surface tempera- 
tures (SSTs) [e.g., Butchart and Scaife, 2001; Austin et al., 

80 2007; Olsen et al., 2007; Li et al., 2008; Eyring et al., 2010a; 
Eyring et al., 2010b; Austin et al., 2010]. However, perform- 
ing numerous sensitivity simulations to separate the different 
chemical processes that control stratospheric changes can be 
more easily done using two-dimensional (2D) models, given 
86 their much smaller computational requirements. 2D mod- 
els have been widely used in international assessments of 
the stratosphere [e.g., WMO, 2003, 2007, 2011], and past 
studies have shown that 2D models can resolve much of the 
large scale stratospheric variability on monthly and longer 
90 time scales, as seen in comparisons with observations and 
3D models [e.g.. Plumb and Mahlman, 1987; Yudin et al., 
2000; Fleming et al., 2007; Newman et al., 2009]. Previous 
2D model studies have investigated the relative roles of the 
long term changes in CO 2 , CH 4 , and N 2 O, focussing on the 
95 stratospheric ozone changes over the next century [e.g., Ran- 
deniya et al., 2002; Chipperfield and Feng, 2003; Portmann 
and Solomon, 2007]. 

In this paper we expand on these previous studies and ex- 
amine in more detail the relative contributions of the tong 
too term changes in atmospheric GHG and ODS loading using 
our recently upgraded Goddard Space Flight Center (GSFC) 
2D coupled chemistry-radiation-dynamics model. We ex- 
ploit the computational speed of the 2D model to perform 
numerous perturbation simulations to investigate the strato- 
105 spheric impacts due to GHG and ODS loading for the 250- 
year time period, 1850-2100. Weexamine the ozone, temper- 
ature, and age of air impacts, and focus on the time periods 
prior to 1950 and the latter half of the 21** century. We also 
use perturbation tests to examine the ozone impacts due to 
no the chemical coupling between CH 4 and chlorine. We then 
investigate the long term time dependence of the photochem- 
ical lifetimes of N 2 O, CFC-1 1 , CFC-12, and CCI4. Here, we 
examine the relative importance of changes in the BDC and 
the photochemical loss rates in controlling these lifetimes. 

115 The recent SPARC Chemistry-Climate Model Validation 
Activity [SPARC CCMVal, 2010] provided a comprehen- 
sive process oriented evaluation of many CCMs. Because 
2D models were not included in this activity, and given 
the recent improvements to our 2D model, we provide in 
120 Appendicies A and B, a detailed description and evalua-os 
tion of our upgraded model, comparing climatological sim- 
ulations with observations of various stratospheric tracers. 
Throughout the paper, we also compare long term simula- 
tions from the 2D model with the Goddard Earth Observing 
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Fig. 1. Time dependent surface boundary conditions for CO 2 , CH4, 
and N 2 O from Hansen and Sato [2004] for 1850-1950 and the 
IPCC GHG .scenario AIB for 1950-2100. The bottom panel shows 
the upper stratospheric Equivalent Effective Stratospheric Chlorine 
(EESC) taken as the global average of Cly-l-OOBry at 50 km. See 
text for details. 


System chemistry-climate model (GEOSCCM) and multi- 
decadal observational data sets to illustrate that the 2D model 
captures the basic processes that drive long term changes 
in stratospheric ozone, temperature, and age of air. The 
good 2D model agreement with the measurements and the 
GEOSCCM then justifies the use of the 2D model for the 
perturbations addressed in this study. 


2 Model Simulations 

For this study, we utilize a series of 2D model experiments in 
which the surface concentrations of only the ODSs or the in- 
dividual GHGs are varied time dependently for 1850-2100, 
while all other source gases are fixed at low (1850) levels. 
In this way, we separate the individual effects of the ODS 
and GHG loading. We compare these with the 2D base- 
line simulation in which all source gases are varied time de- 
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140 pendently, and with the GEOSCCM baseline simulation for 
1950-2100. The GEOSCCM couples the GEOS-4 general 
circulation model with stratospheric chemistry and has been 
applied to various stratospheric problems fe.g., Stolarski et 
al„ 2006; Pawson et al., 2008; Waugh et al., 2009; Oman 
145 et al., 2009; Newman et al., 2009; Li et al., 20091. The 
GEOSCCM uses specified time dependent SSTs and sea-ice 
amounts, and the results presented in this study are com- 
prised of three simulations which utilize somewhat different 
SSTs for the past and future time periods: 1950-2004, 1971- 
,50 2052, and 1996-2100. 

For 1950-2100, the 2D and GEOSCCM simulations use 
surface ODS boundary conditions from scenario Al of 
WMO [2007], and GHG boundary conditions from scenario 
AIB (medium) from the Intergovernmental Panel on Cli- 
155 mate Change (IPCC) Special Report on Emissions Scenarios 
[IPCC, 2000]. The 2D simulations for 1850-1950 use GHG 
surface boundary conditions from Hansen and Sato [2004]. 
For the ODSs , most are zero prior to 1950, except for the fol- 
lowing; CFCI 3 and CF 2 CI 2 are set to zero prior to 1935 and 
,60 1 946, respectively, and are then ramped up slowly to the 1 950 

WMO [2007] values. CCI 4 is ramped up exponentially from 
zero in 1 900 to the 1 950 value of WMO [2007] , approximat- 
ing the time series of Butler etal. [1999]. CHsClandCHsBr 
are set to 440 pptv and 5 pptv respectively in 1 850, and fol- 
165 low the time variation up to 1950 as discussed in Butler et 
al. [1999] and WMO [2003]. The corresponding GHG sur- 
face boundary conditions for 1850-2100 are shown in Figure 
1 , along with the Equivalent Effective Stratospheric Chlorine 
(EESC, bottom panel) representing the time dependent con- 
170 centration of halogen loading from the ODS source gases. 
Here, EESC is taken as the global average of Cly-f 60Bry 
at 50 km. Note that all 2D and GEOSCCM calculations use 
fixed solar flux (no solar cycle variations) and clean strato- 
spheric aerosol conditions specified from WMO [2007]. 


175 3 Ozone 

3.1 Base model-data comparisons 

As a general model evaluation, we first compare the sim- 
ulated vertical profile ozone trends for 1979-1996 with the 
SBUV data for the near-global (60"S-60“N) average (Figure 
180 2, top). These are derived from regression fits to the EESC 
time series (Figure 1) for 1979-2004. At 50 hPa (the lowest 
level of the SBUV data), the observationally-derived trend is 
underestimated in the 2D base simulation (all source gases 
varied time dependently), less so in the GEOSCCM. How- 
,85 ever above 20 km, both base simulations are mostly in rea- 
sonable agreement with the observations. Note that we do 
not show GEOSCCM results in the troposphere as ozone is 
relaxed to an observational climatology in this region. 

The recent past and future ozone changes in the 2D model 
ISO are also mostly similar to the GEOSCCM in the near global 200 
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Fig. 2. (Top) Vertical profiles of the annual and near-global average 
(60"S-60“N) ozone trend (%/decade) for 1996-1979 derived from 
the SBUV data and model simulations. Shown are the base simula- 
tions (all source gases varied time dependently) of the GEOSCCM 
(black dashed lines) and 2D model (black .solid lines), along with 
2D simulations in which only certain source gases are varied as fol- 
lows: ODSs only (blue lines); CO 2 only (red lines); CH 4 only 
(orange lines); N 2 O only (green lines). The trends are derived 
from regression fits to the EESC time series (Figure 1) for 1979- 
2004. (Bottom) As in the top but for the 2095-2005 ozone dif- 
ference (%/decade) using 10-year averages centered on 2095 and 
2005 to reduce the effects of interannual dynamical variability in 
the GEOSCCM. The change due to CH 4 is .shown for 2095-2005 
(orange solid) and 2055-2005 (orange dashed-dotted). 


average (Figure 2) and in the latitude-height variations (Fig- 
ures 3 and 4). This general model agreement is also seen in 
time series of near-global profile ozone (Figure 5), global to- 
tal ozone (Figure 6 ), and tropical total ozone (Figure 7). For 
reference, we include time series of the BUV/SBUV satel- 
lite observations for profile ozone and ground-base data for 
global total ozone (updated from Fioletov et al. [2002]). 

The largest 2D-GEOSCCM differences occur in the 
Antarctic ozone hole region, where the 2D model simulates 
smaller past (negative) and future (positive) ozone changes 
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compared to the GEOSCCM (Figures 3 and 4). Some of this 
is likely due to the 2D model not fully resolving the pro- 
cesses that control polar ozone loss, as these can have large 
zonal asymmetries. Some of this model difference also re- 
205 fleets the known high ozone bias in the high latitude lower 
stratosphere in the GEOSCCM. This bias is most pronounced 
during periods of low chlorine loading so that the chlorine- 
induced changes in the Antarctic spring are too large by 60- 
80% [Pawson et al., 20081. These model differences are re- 
210 fleeted in the near-global averaged vertical profiles below 1 8 
km (Figure 2) and in the total column time series (Figure 6 ) 
in which the GEOSCCM simulates significantly more past 
ozone reduction and future ozone increase. However, the 
generally good agreement between the 2D model base simu- 
215 lations and the observations and GEOSCCM in Figures 2-7 
show that the 2D model captures the basic processes respon- 
sible for long term stratospheric ozone changes. 

3,2 2D perturbation simulations 

The relative roles of ODS, CO 2 , CH 4 , and N 2 O loading in 
220 controlling the recent past and future ozone changes are illus- 
trated by the 2D model perturbation simulations in Figures 
2-7. This includes the well known dominance of ODS load- 
ing in controlling the sharp ozone decline in the lower and 
upper stratosphere globally during ~ 1970-2000. Because of 
225 the strong impact in the upper stratosphere (Figures 3 and 
4), ODS loading also largely controls the tropical total ozone 
time series from ~1970 through the early 21’* century (Fig- 
ure 7). 

CO 2 cooling and subsequent reduction in the ozone loss 
230 rates produce a broad ozone increase of 1 - 2 %/decade in the 
upper stratosphere (Figures 2-4). For the 2005-2095 time 
period, the ozone increases due to increasing CO 2 and de- 
clining ODS emissions are similar in the upper stratosphere 
(1 .5-2%/decade) in Figure 2 (bottom) and Figure 4. By 2100, 

235 CO 2 loading is the dominant impact at 40 km, causing a 
20% increase in ozone from 1850-2100 (Figure 5, middle, 
red curve). 

Previous CCM studies have shown that increased GHG 
loading results in an acceleration of the BDC with related 
240 impacts on ozone and trace gases fe.g., Butchart and Scaife, 
2001; Austin and Li, 2006; Garcia and Randel, 2008; Li et®“ 
al., 2009]. This feature is seen in the 2D C 02 -only sim- 
ulation in Figures 3 and 4, in which lower tropical strato- 
spheric ozone is reduced by 1 - 2 %/decade as ozone-poor air 
24S is advected upwards from the tropical troposphere. There 
is a compensating downward advection of ozone-rich air®“ 
in the extratropics at 10-15 km which is strongest in the 
NH. This hemispheric asymmetry is consistent with previ- 
ous GEOSCCM simulations of the climate change impacts 
250 on ozone [Olsen et al., 2007; Li et al., 2009; Waugh et al., 
2009]. These BDC-driven ozone changes are reflected in theaes 
global average in Figure 2, with an ozone decrease (increase) 
of ~ 0.5%/decade centered near 20 km (15 km) in the 2D 
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Fig. 3. Annually averaged ozone trend for 1996-1979 derived from 
the GEOSCCM and 2D model simulations. The top panels are from 
the ba.se simulations in which all source gases are varied time de- 
pendently. The middle and bottom panels show 2D simulations in 
which only the ODSs or GHGs are varied time dependently as indi- 
cated. The trends are derived from regression fits to the EESC time 
series (Figure 1) for 1979-2004. The GEOSCCM trends are not 
computed in the troposphere as ozone is relaxed to a climatology in 
this region. The contour intervals are ±2 %/decade and include the 
±.5 and ±1 %/decade contours. 


C 02 -only simulation. The COa-induced ozone decrease is 
dominant by 2100 in the near global average time series at 
22 km (Figure 5 , bottom, red curve) , with a decrease of 5 .5% 
from 1850-2100. At this level, N 2 O and CH 4 loading have 
negligible impacts, so that the net result of ODS and CO 2 
changes is an ozone increase from 2000-2030 and a decrease 
from 2030-2100 in the base simulations. 

The net impact of CO 2 loading on total ozone is an in- 
crease of 12.5 DU (4.2%) from 1850-2100 in the global av- 
erage (Figure 6 , red curve). CO 2 loading also increases total 
ozone at midlatitudes of both hemispheres (not shown), and 
as with profile ozone, the total column increase is more pro- 
nounced in the NH compared to the SH midlatitudes owing 
to the larger enhancement of the BDC north of the equator. 
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Fig. 4. As in Figure 3, except for the 2095-2005 ozone difference 
using 10-year averages centered on 2095 and 2005 to reduce the ef- 
fects of interannual dynamical variability in the GEOSCCM. The 
contour intervals are ±2 %/decade and include the ±.5 and ±1 
%/decade contours. 


In the tropics, the enhancement of the BDC advecting ozone- 
poor air from the troposphere counteracts the ozone increase 
Z70 in the upper stratosphere caused by the CO2 cooling. The net 
effect is a decrease in tropical total column ozone through- 
out the 21’* century in Figure 7 (red curve), consistent with 
previous CCM results [Li et al., 2009; Waugh et al., 2009; 
Eyring et al., 2010b]. 

275 In Figures 2-5, N2O loading and the subsequent increase 
in stratospheric NOy lead to a decrease in stratospheric 
ozone, with a maximum decline of —0.5 to —0.6%/decade 
near 35 km in the global average in both the past and future. 290 
These magnitudes are slightly less than obtained by Port- 
280 mann and Solomon [2007] who used the IPCC A2 GHG sce- 
nario which has larger N2O increases compared to the AIB 
scenario used here. From 1850-2100, N2O loading results in 
a total ozone decrease of 8 DU (—2.7%) in the global average zas 
(Figure 6) and 4 DU (— 1 .5%) in the tropics (Figure 7). This 
285 effect of N2O loading taken in isolation is larger than would 
be when taking into account the effects of increased CO2 
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Fig. 5. Near global (60“S-60'’N) annually averaged ozone time se- 
ries for 1860-2100, relative to 1860, at 22, 40, and 60 km. Values 
are in ppmv change (left axes) and % change (right axes). Shown 
are the base simulations (all source gases varied time dependently) 
of the GEOSCCM (black dotted lines) and 2D model (black solid 
lines), along with 2D simulations in which only certain source ga.ses 
are varied as indicated. The GEOSCCM time series has been ad- 
justed to match the 2D base simulation for 1960. Also shown are the 
BUV/SBUV satellite observations for 1970-2009 (excluding 1973- 
1978) at 22 and 40 km (”+” symbols). To emphasize the model-data 
comparison after 1 970, the data have been adjusted so that the 1970- 
1 972 average matches that of the base simulations. 


cooling on future NOy and ozone changes, as has been dis- 
cussed previously [Rosenfleld and Douglass, 1998; Daniel et 
al., 2010]. We will discuss this in more detail in section 3.5. 

Atmospheric CH4 impacts ozone via three mechanisms: 
1) increases in CH4 increase the amount of H2O in the strato- 
sphere and mesosphere which in turn reduces ozone by en- 
hancing the HOx -ozone loss cycles; 2) increases in CII4 
lead to increased ozone throughout the stratosphere by con- 
verting active chlorine to the reservoir HCl via the reaction 
CH4 + CI-7HCl-fCH:i; and 3) increases in CH 4 lead to 
increased ozone in the troposphere and very lower strato- 
sphere due NOx-induced ozone production [e.g, Brasseur 
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Fig. 6. Globa! and annual averaged total ozone time series relative 
to 1860 values. Shown are the base simulation (all source gases 
varied time dependently) from the 2D model (black solid line) and 
GEOSCCM (black dotted line). The GEOSCCM time series has 
been adjusted to match the 2D base simulation for 1960. Also 
shown are 2D simulations in which only certain source ga.ses are 
varied time dependently as indicated, with the other source gases 
fixed at 1850 levels. The ”+” symbols represent ground-based data 
updated from Fioletov et al. [2002]. To emphasize the model-data 
comparison after 1970, the data have been adjusted so that the 1964- 
1970 average matches that of the base simulations. Also shown 
are the combined effects of CH4 and ODSs loading: the 2D-CH4 
only(a) curve (orange dashed-dotted) is the difference between a 
simulation with CH4 and the ODSs varied time dependently and 
that with only the ODSs varied time dependently, showing the effect 
of CH4 in the presence of time dependent ODS loading; the 2D- 
ODSs only(b) curve (blue dashed-dotted) is the difference between 
a simulation with CH4 and the ODSs varied time dependently and 310 
that with only CH4 varied time dependently, showing the effect of 
ODSs in the presence of time dependent CH4 (CO2 and N2O are 
fixed at 1850 levels in these simulations). Values are in Dobson Unit 
(DU) change (left axes) and % change (right axes). 


and Solomon, 1986]. This latter mechanism is strongly de- 
300 pendent on the amount of tropospheric NOx- To more prop- 
erly account for this in the 2D model, we constrain the model 


Fig. 7. As in Figure 6 except for tropical total ozone (10°S-10°N 
average), and without the 2D-CH4 only(a) and 2D-ODSs only(b) 
curves. For visual clarity, the observations are not shown as these 
are dominated by the large quasi-biennial oscillation, a feature not 
included in the simulations. Values are in Dobson Unit (DU) change 
(left axes) and % change (right axes). 

troposphere using output from the Global Modeling Initia- 
tive’s (GMI) combined stratosphere-troposphere chemistry 
and transport model (GMI Combo CTM) [Strahan et al., 
305 2007; Duncan et al,, 20071. This includes specifying the 
lower tropospheric (below 5 km) NOy boundary condition 
to account for surface emission and loss processes, NOx pro- 
duction from lightning, and washout of HNO3 , all for present 
day conditions. The resulting 2D tropospheric NOx and 
NOy distributions compare favorably with the GMI model. 
The 2D model utilized in this way simulates a mid-upper tro- 
pospheric ozone response to a steady state CH 4 perturbation 
(0.5 ppmv) similar to that obtained in the GMI model for 
present day conditions. 

Figures 2-5 show that the combined effect of the three 
CH 4 mechanisms outlined above yield a significant ozone 
increase (1-1 .5 %/decade) in the troposphere and very lower 
stratosphere, a mostly weak positive response at 20-40 km, 
and a negative ozone response above 40 km with decreases 
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320 of — 1%/decade and larger above 60 km. Figure 5 (top panel) 
shows that CH 4 loading is dominant in controlling the ozone 
time dependence at 60 km, due to the subsequent increase in 
H 2 O and HOx-ozone loss. This mechanism also results in 
significant CH 4 -induced ozone loss at 60 km prior to 1960 
325 so that one-third of the total decrease in ozone during 1 850- 
2050 has occurred by 1960. For total ozone, CH 4 loading 
results in increases from 1850-2050 of 7.5DU (2.5%) in the 
global average (Figure 6 ) and 4.5 DU (1.7%) in the tropics 
(Figure 7). 

330 The future ozone response to CH 4 follows the surface 
boundary condition, which is significantly larger in 2055 
than in 2100 in the AIB scenario (Figure 1). This results 
in a relatively small change from 2005-2095 seen in the 
latitude-height section in Figure 4 (< ±0 .5%/decade every- 
335 where). We show vertical profile results for both 2005- 
2055 and 2005-2095 in Figure 2, illustrating the much dif- 
ferent magnitudes of the response for the different time peri- 
ods. These responses are also qualitatively similar but have 
smaller magnitudes compared to Portmann and Solomon 
340 [2007] who used IPCC scenario A2 which has larger GHG 

loading. We note that because the 2D model tropospheric 
NOx and NOy are fixed to the present day GMI simulation, 
the tropospheric ozone response to the time dependent CH 4 
perturbation shown in Figures 2-7 does not account for long 
345 term changes in tropospheric NOx emissions. These have 
undergone significant past increases which would lead to 
more CH 4 -induced tropospheric ozone production than we 
show in Figures 2-7. 

By 2100, ODS, CO 2 , CH 4 , and N 2 O loading all play im- 
350 portant roles in controlling global total ozone. CO 2 has the 
largest impact, leading to a 4% increase from 1860-2100, 
with CH 4 loading causing a 2% increase, and N 2 O and ODS 
loading each contributing a 2-2.5% decrease from 1860-2100 
(Figure 6 ). The net result is a 1 .7% (5DU) increase from 
355 1 860-2 1 00 in the 2D base simulation . In the tropics , baseline 

total ozone increases during 2000-2050. This is followed 
by a ~ 1 % decrease over the second half of the century, as 
the combination of increasing CO 2 and N 2 O and decreas- 
ing CH 4 more than offset the effect of reduced ODS load- 
360 ing. This general time variation of 21 century tropical total 
ozone is consistent with the GEOSCCM in Figure 7 and other 
CCM results reported recently [Eyring et al., 2010a, 2010b; 
Austin et al., 2010]. 

33 1860-1960 Ozone 

365 As seen in Figures 5-7, the ozone changes prior to 1960 
are relatively small, but not insignificant, compared to post- 
1960. To examine this more closely. Figure 8 shows a close- sao 
up of the 1860-1960 ozone time series for the near-global 
average at 40 km and the global average total column. 

370 From 1860-1960, the total column impacts due to CO 2 
and N 2 O are approximately equal and opposite, being 
-bO.5% and —0.5%, respectively. The CH 4 impact is slightly sss 



1.0 


Global total column 


V 

O' 

c 


3 

o 


-2 

-3 


- 2 D-bose 

- 2 D- 0 DSS only 
-2D-CO2 only 

■ 20-CH4 only 
2D-N2O only 



<u 

O’ 

c 




0.5 


- 1.0 

1860 1880 1900 1920 1940 1960 


Fig. 8. Ozone time series for 1860-1960, relative to 1860 values, 
from 2D model simulations for the 60°S-60°N average at 40 km 
(top), and the 90°S-90"N average total column (bottom). Shown 
are the baste simulation (all source gases varied time dependently) 
and simulations in which only certain source gases are varied time 
dependently as indicated, with the other source gases fixed at 1860 
levels. Values are in Dobson Unit (DU) change (left axes) and % 
change (right axes). 


larger by 1960 (-1-0.75%). The ODS impact, which is due 
mainly to CCI 4 emissions, causes a 1 % depletion in the total 
column by 1960, with the vast majority of this decrease oc- 
curing after 1920. At 40 km, the impact of CO 2 cooling leads 
to a 2% ozone increase from 1860-1960, which is approxi- 
mately equal and opposite to the ozone depletion caused by 
ODSs. 

Ozone in the resulting base simulation (black curves) 
reaches broad pre-modem maxima during 1920-1940, with 
increases from 1860 of 0.8 DU (0.3%) in the total column 
and 0.06 ppmv (0.8%) at 40 km. This preceeds the decline in 
ozone driven mainly by ODS loading, which becomes much 
more rapid after ~ 1970 (Figures 5 and 6 ). 
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3.4 CH4 sensitivity experiments 

As the changes in CH4 impact ozone via the different mecha- 
nisms outlined above, we now examine these effects in more 
detail. Figure 9 (top panel) shows the global/annual aver- 
aged 2D model steady state ozone change due to a 0,5 ppmv 
CH4 perturbation for year 2000 conditions. This is expressed 
in DU/km to emphasize the contribution to the total column 
change. CH4 is lo.st via reaction with 0(^D), OH, and Cl, 
and the black curve shows the standard case in which all three 
reactions use the perturbed CH4. This shows that the CH4 
impact on ozone is positive everywhere below 46 km, with 
only a very small negative change above 46 km, and a total 
column response of +2.7DU. 

To qualitatively separate the effects of the different mech- 
anisms, the green curve is a simulation in which perturbed 
CH4 is used for the reactions with 0(*D) and OH, with un- 
perturbed CH4 used for the reaction CH4 -l-Cl. In this way, 
the conversion of active chlorine (Cl) to the reservoir (HCl) 
is not impacted by the additional CH4 . Therefore, ozone only 
responds to the enhanced NOx-induced ozone production in 
the troposphere yielding ozone increases in this region, and 
the enhanced H2O and HOx-ozone loss cycles above ~35 
km and in the very lower stratosphere, which yield ozone de- 
creases. The enhanced H2O in this case also enhances polar 
stratospheric clouds and the heterogeneous chemical ozone 
loss in the SH polar region, which results in a negative re- 
sponse at 17-28 km in the global average in Figure 9. 

The red curve in Figure 9 shows the opposite case of the 
green curve, i.e., here, perturbed CH4 is used for the reac- 
tion CH4 -I- Cl, with unperturbed CH4 used for the reactions 
with 0(^D) and OH. This significantly reduces the amount 
of H2O and HOx-ozone loss, as well as the NOx -ozone pro- 
duction in the troposphere generated by the CH4 perturba- 
tion. This better isolates the impact of CH4 in controlling 
the chlorine partitioning and subsequent chlorine-catalyzed 
ozone loss. The resulting ozone response is positive through- 
out the stratosphere, with the largest impacts occurring in 
the regions where chlorine-catalyzed ozone destraction is 
largest, i.e., the upper stratosphere globally, and in the lower 
stratosphere corresponding to the ozone hole region. Com- 
paring the green and red curves shows that the CH4 + Cl re- 
action has the largest impact on ozone at ~ 15-45 km, with‘>« 
the NOx -ozone production mechanism being dominant be- 
low 15 km. 

The bottom panel in Figure 9 shows the same cases run for 
2100 conditions, which has greatly reduced chlorine loading 
compared with present day conditions. The effect of the fu-44s 
ture reduced chlorine loading is evident, as the case in which 
the perturbed CH4 is used only for the CH4-I-CI reaction 
(red curve) exhibits a much weaker ozone response at all lev- 
els compared with the simulations for present day conditions. 
For the total column, this simulation accounts for only 28% 450 
of the full CH 4 response (+1.8DU) in 2100, compared to 
70% of the full response in 2000. In 2100, the full ozone 


Ozone change due to CH4 perturbation 

Global/Annual average 




+ 1.8DU 

-H.30U 

+0.50U 


.12 


60 

50 

40 

;30 

;20 

10 

0 

.16 


— All CH4 reactions perturbed 

— only CH 4 +OH and CH 4 + 0 'D perturbed 

— only CH 4 +CI perturbed 


Fig. 9. Vertical profiles of the annual and global averaged steady 
state ozone change (DU/km) due to a 0.5 ppmv CH4 perturbation 
for year 2000 (top) and year 2100 (bottom). The black curves show 
simulations with all CH4 reactions (reactions with 0('D), OH, and 
Cl) using perturbed CH4 . The green curves show simulations us- 
ing perturbed CH4 for reactions with 0(*D) and OH, with unper- 
turbed CH4 used for the reaction CH4 + Cl. The red curves show 
simulations using perturbed CH4 for the reaction CH4+CI with 
unperturbed CH4 used for reactions with 0(*D) and OH. The total 
column responses for each case are listed in the upper right corner 
of each panel. 


response to the CH4 perturbation (black curve) is dominated 
below ~I8 km by the enhanced NOx-ozone production, and 
above ~45 km by the enhanced HOx-ozone loss. 

For time dependent total ozone, the combined effect of 
CH4 and Cly loading is illustrated in Figure 6. The solid or- 
ange curve (2D-CH4 only) uses fixed 1 850 (very low) chlo- 
rine loading, and is therefore similar to the green curves in 
Figure 9 (unperturbed CH4 + CI) since the impact of the 
CH4-I-CI reaction is minimal in this case. To illustrate the 
full effect of CH4 in the presence of time dependent ODS 
loading (analogous to the black curves in Figure 9), the 2D- 
CH4 only (a) curve (orange dashed-dotted in Figure 6) is the 
difference between a simulation with CH4 and the ODSs var- 
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ied time dependently and that with only the ODSs varied time 
dependently (CO 2 and N 2 O are fixed at 1850 levels in both 
455 simulations). Compared to the simulation using fixed 1850 
Cly (solid orange line), the full effect of CH 4 loading results 
in significantly more global total ozone when the ODS load- 
ing is highest, i.e., ~1985 through the first half of the 21'’* 
century. 

460 This effect can also be seen in the ODS-only simulation, 
i.e., the effect of ODS loading in the presence of time depen- 
dent CH 4 (the 2D-ODS only (b) curve (blue dashed-dotted) 
in Figure 6 ). This is the difference between the simulation 
with CH 4 and the ODSs varied time dependently and that 
465 with only CH 4 varied time dependently. Comparing the two 
blue curves illustrates how the impact of the ODS loading is 
mitigated by the time dependent CH 4 changes. This effect is 
largest when the ODS loading is largest, i.e., ~2000, when 
the global total ozone depletion due to ODS loading is 3 DU 
470 ( 1 %) less in the presence of the time dependent CH 4 changes 

compared to that with fixed 1850 CH 4 levels. 

3,5 Quantification of the Ozone Impacts due to GHGs 

Given the importance of increasing GHG emissions on fu- 
ture ozone levels, we examine these relative impacts in more 
475 quantitative detail. To do this, we apply the concept of the 
Ozone Depletion Potential (ODP) to GHG emissions, as re- 
cently done for N 2 O and the induced NOx-ozone destruc- 
tion [Ravishankara et al., 2009; Daniel et al., 2010]. Tra- 
ditionally, the ODP metric has been used to quantify the 
460 change in global ozone per unit mass emission of a spe- 
cific chlorine-containing compound relative to the change 
in global ozone per unit mass emission of CFC- 1 1 (CFCl,-j) 
[Wuebbles, 1983; Fisher et al., 1990; Solomon et al., 1992]. 
Because the ODP is defined as a ratio to the ozone loss due tosio 
485 CFC-11, many uncertainties in the ozone-loss computation 
cancel. For this reason ODPs of chlorocarbon compounds 
are generally less sensitive to photochemical modelling un- 
certainties than are absolute ozone loss calculations. 

Application of the ODP concept to a non-chlorine con-s's 
490 taining compound does not benefit in the same way from the 
cancellation of uncertainties in the ratio of ozone loss to that 
of CFC-1 1, e.g. N 2 O does not cause an ozone hole. More 
generally, the chemical mechanisms that impact ozone are 
different for GHGs, and tend to occur in different regions 520 
495 of the atmosphere, compared to those of chlorine containing 
compounds. 

While being cautious about the interpretation of the ODP 
concept as applied to non-chlorine containing compounds, 
we use the ODP metric to examine the relative impacts on 525 
soo ozone of N 2 O, CH 4 , CO 2 , and CFC-11 for 1850, 1950, 
2000, and 2100 steady-state atmospheric conditions. These 
calculations are meant as a guide for evaluating the relative 
importance of GHG loading on past and future ozone levels. 

We use the standard method to calculate ODPs, i.e. ,530 
505 change each species by an amount that leads to a 1 % change 


Table 1. Steady state 2D model calculated annual/global average 
flux (kg/sec) required to produce a 1 % change in annually averaged 
global total ozone for the compounds and years listed. The ozone 
change is negative for CFC-1 1 and N2O, and positive for CH4 and 
CO 2 . 


1850 

1950 

2000 

2100 

CFC-11 1.52 

1.65 

2.02 

1.83 

N 2 O 61.9 

68.5 

107 

104 

CH 4 4980 

4960 

3130 

5960 

CO 3 241 

269 

519 

881 

Table 2. Steady state 2D model calculated ozone depletion potential 
(ODP) relative to CFC-1 1 for N20,CH4, and CO 2 for the years 
indicated, based on the values listed in Table 1 . 

1850 

1950 

2000 

2100 

N 2 O -1-0.025 

-f 0.024 

+0.019 

+0.018 

CH 4 -0.00031 

-0.00033 

-0.00065 

-0.00031 

CO 2 -0.0063 

-0.0062 

-0.0039 

-0.0021 


in annually-averaged global total ozone. Table 1 shows the 
resulting annual and global average flux of each compound 
for the years indicated. The ozone change is negative for 
CFC-1 1 and N 2 O, and positive for CH 4 and CO 2 . Note that 
the flux changes for N 2 O, CH 4 and CO 2 are much larger, 
by one to three orders of magnitude, than those for CFC- 
1 1 . Table 2 shows these numbers converted to the standard 
ODP, i.e., the ratio of the flux change needed to cause a 1% 
ozone change for each compound to the flux change for CFC- 
11. Our N 2 O ODP values for 1950 and 2000 (-b0.024 and 
-FO.019, respectively), are similar to those reported in Ravis- 
hankara et al. [2009]. These results also illustrate that while 
N 2 O has a positive ODP, CH 4 and CO 2 have negative ODPs 
relative to CFC-1 1 . However because the fluxes and ODPs 
vary by two orders of magnitude, and the background con- 
centrations also undergo large changes, it is difficult to eval- 
uate these values on an absolute basis or relative to one an- 
other. 

Tables 3 and 4 illustrate another way to look at the problem 
of comparing the impact of GHGs on ozone. Table 3 shows 
the sensitivity of total column ozone to a specified change in 
the concentration of each GHG. These are expressed in per- 
cent change in ozone for a 1 ppbv or 1 ppmv change in the 
surface concentration of the GHGs. The impact of each can 
then be deduced by multiplying by the actual change in the 
concentration of the species over a given time period. This is 
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535 


540 


545 


550 


555 


Table 3. Steady state 2D model calculated percentage change in 
annual/global average total ozone per unit mixing ratio change for 
the compounds and years listed. The mixing ratios are in parts per 
billion by volume (ppbv) for N 2 O, and parts per million by volume 
(ppmv) for CH4 and C02- Shown in parentheses are the percent- seo 
age changes of each compound, relative to the background levels, 
required to produce a 1% ozone change. 



1850 

1950 

2000 

2100 

N 20 

-0.031 (12) 

-0.028(13) 

-0.020(16) 

-0.021 (13) ® 

CH4 

-H.8 (70) 

+ 1.7 (49) 

+2,6 (22) 

+ 1.4 (36) 

C 02 

-1-0.016(23) 

+0.014(23) 

+0.0077 (34) +0.0041 (33) 


570 


Table 4. Percentage change in annual/global average total ozone 
for the compounds and time intervals listed, based on the sensitivity 
factors in Table 3 and the change in the surface boundary condition 
in the AIB scenario (Figure 1). The CH 4 -induced ozone change is 
negative for 2080-2100 due to the decrease in the methane boundary 
condition during this time. 



1940-1960 

1980-2000 

2080-2100 

N 2 O 

-0.12 

-0.30 

-0.15 

CH 4 

+0.3! 

+0.49 

-0.30 

CO 2 

+0.07 

+0.25 

+0.34 


similar to the method used in Stolarski et al. [20101 for tem- 
perature changes. The resulting percentage ozone changes 
for several 20-year periods are listed in Table 4. Note that the 
CH 4 -induced ozone change is negative for 2080-2100 since 5 „ 
the methane boundary condition decreases during this time 
period (Figure I). These calculations illustrate that CII4 had 
the largest GHG impact on ozone during 1980-2000, owing 
to the large effect of high Cly loading, with the N 2 O and 
CO 2 impacts being roughly equal and opposite. The CH4 5,5 
impact was also substantially larger than either N 2 O or CO 2 
during 1940-1960, reflecting the large relative changes in the 
CH 4 boundary condition during this time. However at the 
end of the 21 ** century, the CO 2 impact is projected to be 
the largest, and will be more than twice the magnitude of thCsoo 
N 2 O impact. 

Table 3 also shows the sensitivities as percent changes to 
the background levels of each GHG, required to produce a 
1% ozone change (listed in parentheses). For N 2 O, a ~15% 
change in the background is required to cause a 1 % change in 
ozone, compared to a 4^25-35% change for CO 2 , and ~ 20 - 
70% change for CH 4 . ms 

The time dependence of the values in Tables I -3 partly re- 
flects the fact that as the background concentrations increase, 
more flux of a substance is required to produce a 1 % ozone 


change. However, this time dependence also reveals interest- 
ing aspects of the interactions of the perturbations. For CFC- 
1 1, increasing levels of N2O (NOy) and CH 4 convert active 
chlorine to the reservoirs CIONO 2 and HCl, respectively, 
thereby reducing the efficiency of ClOx-ozone destruction. 

For N 2 O, the time dependence is partly due to the chang- 
ing background Cly, i.e., with higher Cly more NO^ is tied 
up in CIONO 2 , thereby reducing the efficiency of the NOx- 
ozone loss. Other factors include: 1) the CO 2 cooling of 
the stratosphere which results in greater chemical destruc- 
tion of NOy [Rosenfleld and Douglass, 1998], and therefore 
reduced NOx — Orj loss; and 2) a general decrease in 0(^D), 
by as much as 10%, from ~I970-2100 throughout most of 
the stratosphere which results in less NOx production via the 
reaction N 20 -l- 0 (*D). As a result of these factors, our N 2 O 
ODP in 2100 is a bit smaller than in 2000. The N 2 O ODPs 
for 1850 and 1950 are very similar, as are the CH4 and CO2 
ODPs, reflecting the small GHG and EESC-induced changes 
in the stratosphere prior to 1950 (Figures 5 and 6 ). 

The ozone impact of CH4 is strongly dependent on the at- 
mospheric chlorine loading, via the reaction CH 4 -f Cl which 
affects the partitioning of Cly (e.g.. Figures 6 and 9). As 
a result, significantly less CH4 flux is required (Table 1) to 
get a 1% ozone increase in 2000 (large Cly loading) com- 
pared to 1850, 1950, and 2100 (small Cly loading). Sim- 
ilarly, the CH4 perturbation as a percentage of the back- 
ground (Table 3) is significantly smaller in 2000 compared 
to the other years. The resulting CH 4 ODP is nearly twice 
the magnitude in 2000 compared to 1950, 1950, or 2100. As 
discussed earlier, the CH4 results presented here do not ac- 
count for changes in tropospheric NOx emissions (the model 
uses fixed present day NOy specified from the GMI sim- 
ulations), These likely undergo significant past and future 
changes [IPCC, 2007], with less (more) NOx leading to less 
(more) tropospheric ozone production and a smaller (larger) 
negative CH4 ODP than we report in Table 2. 

For CO2, the strong time dependence in Tables 1-3 is 
mostly due to the changing background concentration, which 
increases by a factor of 2.5 from 1 850-2100 (Figure 1 ). How- 
ever, the percentage change relative to the background in Ta- 
ble 3 is somewhat larger in 2000 and 2100 than in 1850 and 
1950. This is partly due to the fact that CO 2 cooling and the 
corresponding increase in upper stratospheric ozone become 
less efficient at higher CO 2 levels, i.e., a saturation effect 
[e.g., Ramaswamy et al., 2001a]. There are other higher or- 
der effects that influence this CO 2 time dependence, includ- 
ing the future changing PSC concentrations - due to strato- 
spheric cooling and increasing water vapor - coupled with 
the decreasing halogen loading. 

4 Temperature 

Previous studies have investigated stratospheric temperature 
trends in the recent past and future, including the relative 
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Fig. 10. Global and annually averaged temperature time series for 
1850-2100 at 20, 40, and 60 km from the base simulations (all 
source gases varied time dependently) of the GEOSCCM (black 
dotted lines) and 2D model (black solid lines). Also shown are 2D 
simulations in which only certain source gases are varied time de- 
pendently as indicated. The right hand axes show the changes rel- 
ative to 1 850. The observations are from the NCEP reanalysis and 
reanalysis-2 data at 20 km for 1958-2009 (blue ”-h”), and the new 
NASA Modem Era Retrospective-analysis for Research and Appli- 
cations (MERRA) meteorological analyses for 1979-2009 at 40 and 
60 km (purple ”+”). To account for the systematic differences be- 
tween the models and data, we have added the following offsets 
to the models at 20, 40, and 60 km, respectively: 2D model: OK, 
-3.4K, -1.5K; GEOSCCM: -2.7K, -H .3K, -2K. 


contributions of ODS and GHG loading [e.g., Ramaswamy 
et al., 2001b; Shine et al., 2003; Shepherd and Jonsson, 

610 2008; Stolarski et al., 20101. Here we briefly summarize the 
model simulated temperature changes, focussing on the rela- 
tive contributions of CO 2 , CH 4 , and N 2 O. 

Figure 10 shows the globally and annually averaged tem- 
perature time series from the base GEOSCCM and 2D 
615 model simulations as indicated. As an observational ref- 
erence, we also show the NCEP reanalysis and reanalysis - ks 
2 data at 20 km for 1958-2009, and the new NASA Mod- 


Fig. 11. (Top) Vertical profiles of the annual and global averaged 
temperature change (K/decade) for 20(X)-1960 (top) and 2095-2005 
(bottom) from the 2D model simulations as indicated. The change 
due to CH4 is shown for 2095-2005 (orange solid) and 2055-2005 
(orange dashed-dotted) in the bottom panel. 

ern Era Retrospective-analysis for Research and Applica- 
tions (MERRA) meteorological analyses (see the wesite: 
620 http://gmao,gsfc.nasa.gov/research/merray) for 1979-2009 at 
40 and 60 km. To match the early years of the data and to ac- 
count for the systematic differences between the models, the 
model curves have been offset as listed in the Figure 10 cap- 
tion. Although there is interannual variability in the data, the 
625 base models simulate the general rate of stratospheric cooling 
seen in the observations. Also, the 2D and GEOSCCM base 
simulations have similar rates of cooling throughout 1950- 
2100, consistent with the fact that both models utilize the 
same IR radiative transfer schemes. Further model tempera- 
630 ture trend comparisons are discussed in Appendix B (Figures 
B2 and B3). 

Following the time change in the surface boundary con- 
ditions (Figure 1), the 2D simulated temperature changes in 
Figure 10 are significantly larger after ~ 1970 compared with 
the 1850-1970 time period. Prior to 1950, the temperature 
changes in the stratosphere are due almost entirely to the 
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CO 2 loading, with negligible impacts due to ODS and N 2 O 
loading. The temperature changes for 1850-1950 from the 
2D base run are mostly linear at 20, 40, and 60 km, andeso 
640 are — 0.05K/decade, — 0.14K/decade, and -0.17K/decade, 
respectively. These magnitudes are generally consistent with 
the GISS modelE results reported in Hansen et al. [2007]. 

As discussed in previous studies, ODS loading and the cor- 
responding reduction in ozone heating have a large impact on 
645 the temperature changes in the recent past, maximizing in the 
lower and upper stratosphere as seen in the vertical profiles 
in Figure 11 (top). By 2100, CO 2 cooling dominates the 
total temperature change throughout the stratosphere, and is 
significantly larger than the warming caused by the reduced™" 
650 ODS loading (Figure 1 1 , bottom). 

The temperature impact of N 2 O caused by the associated 
NOx-ozone depletion maximizes at 30-40 km but is rela- 
tively small (-0.05 K/decade in Figure 11), corresponding to 
a cooling of -0.7K from 1850-2100 at 40 km (green curve in™® 
656 Figure 10, middle). 

Increases in CH 4 Induce temperature changes via several 
mechanisms. Below ~1 8 km, there is a slight warming (Fig- 
ure 1 1 ) caused by the combination of the direct infrared ra- 
diative effect of CH 4 and the indirect radiative effect of the 
660 resulting ozone increases in the troposphere and lower strato- 
sphere [e.g, Portraann and Solomon, 2007]. In the strato- 
sphere, increases in CH 4 lead to warming via the increases 
in ozone caused by the reduction of active chlorine by the 
the reaction CH 4 -F Cl — + HCl + CH3 . Increases in CH4 also^'° 
665 lead to cooling in the stratosphere and mesosphere via the 
increases in water vapor. In the global average, the combina- 
tion of these processes results in a net cooling above ~20 km 
in Figures 10-1 1. 

In the lower mesosphere, the cooling effect of the CH 4 -”° 
670 induced H 2 O increases becomes large and is comparable to 
the CO 2 cooling during 1850-2000 at 60 km (Figure 10, top). 

By 2055, the CH 4 -induced cooling is roughly half that due 
to the CO 2 changes at 60 km in Figure 10. We note also 
that the future CH 4 -induced temperature changes in Figure 
675 1 1 (bottom) are significantly smaller by 2095 compared to 

2055, following the decrease in the CH 4 surface boundary 
condition during the latter half of the 2D* century in the AIB 
scenario (Figure 1). 


5 Age of Air 

660 Previous studies have shown that GHG loading and ozone 735 
depletion impact the stratospheric circulation and age of 
air [e.g., Butchart and Scaife, 2001; Austin and Li, 2006; 
Butchart et al., 2006; Kodama et al., 2007; Garcia and Ran- 
del, 2008; Li et al., 2008; Oman et al., 2009]. Although ob- 

685 servations show little or no age trend over the past 30 years 740 
in the NH above 24 km [Engel et al., 2009], simulation of a 
general decrease in age through the 2 1 ** century appears to 


be robust among coupled chemistry-climate models [SPARC 
CCMVal,2010], 

Here we briefly examine how the individual source gas 
perturbations used in this study impact the model age of air 
simulations. This is summarized in Figure 12 (top panel), 
which shows time series of global and annually averaged age 
of air at 25 km from the various model simulations as in- 
dicated. These results are representative of altitudes above 
~20 km. There is significant interannual variability in the 
GEOSCCM time series which is comprised of three simula- 
tions covering the periods 1950-2004, 1971-2052, and 1996- 
2100. However, the base 2D and GEOSCCM simulations 
have similar rates of decrease in mean age over the 1960- 
2100 time period. 

As with ozone and temperature, the 2D-simulated changes 
in mean age are relatively small prior to 1950, with CO 2 
loading accounting for about half of the total decrease of 0.1 
years from 1860-1950. The base simulated age decreases 
much more rapidly after ~1970 when the impacts of CO 2 
and ODS loading become large and act in the same direc- 
tion. The total age decrease is 0.3 years from 1980-2005, 
with 65% of this change due to ODS loading, 25% due to 
CO 2 , and 10% caused by N 2 O and CH 4 changes. This dom- 
inance of the ODS loading is consistent with the GEOSCCM 
analysis of Oman et al., [2009]. Figure 12 (top) also shows 
a 2D-ODS only simulation with all heterogeneous chemical 
processes turned off, i.e., no ozone hole is simulated (blue 
dashed-dotted line). This simulation shows only small mean 
age change due to gas-phase chlorine-ozone destmction, and 
suggests that most of the ODS effect on the age of air change 
is due to the ozone hole. Separation of these ODS impacts on 
mean age is currently being investigated in the GEOSCCM 
(L. Oman, personal communication). 

In the future, the base simulation age decreases less 
rapidly as the effects of the reduced ODS loading partially 
offset the CO 2 impacts, with the latter becoming the dom- 
inant mechanism after about 2025. By 2100, CO 2 loading 
accounts for ^75% of the total age change (~1 year) from 
1860. N 2 O, CH 4 , and the remaining effects of ODS loading 
have secondary impacts by 2100 , with mean age decreases 
of 0.12, 0.05, and 0.08 years (12%, 5%, and 8 % of the total) 
respectively, from 1860-2100. 

Previous studies have attributed the BDC acceleration to 
increased stratospheric wave driving resulting from changes 
in the zonal mean winds. These are ultimately due to changes 
in the temperature distribution induced by SST changes, 
GHG loading, and polar ozone depletion [e.g., Olsen et al., 
2007; Kodama et al., 2007; Garcia and Randel, 2008; Oman 
et al., 2009]. The 2D model formulation uses parameter- 
izations to account for planetary and gravity wave effects 
(section Al), and the latitude-height patterns of the long 
term changes in zonal mean temperature, zonal wind, wave- 
induced acceleration, lower stratospheric tropical upwelling, 
and age of air are generally similar to those in WACCM3 and 
GEOSCCM [Garcia and Randel, 2008; Oman et al., 2009]. 
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Fig. 12. Global and annually averaged age of air at 25 km. (Top) 
Time series for 1860-2100 from the base simulations (all source 
gases varied time dependently) of the GEOSCCM (black dotted 
lines) and 2D model (black solid line.s). Also shown are 2D simula- 
tions in which only certain source gases are varied time dependently 
as indicated. The blue dashed-dotted line is a 2D-ODS only simu- 
lation with all heterogeneous chemical reactions turned off so that 
no ozone hole is simulated. The right hand axis shows the change 
relative to 1860. The GEOSCCM time series is comprised of three 
simulations, which use somewhat different SSTs, for the time peri- 
ods 1950-2004, 1971-2052, and 1996-2100 (the first 8-10 years of 
each simulation have been removed to allow for spin-up). For visual 
clarity and to account for a systematic offset between the models, 
we have added .3 years to the GEOSCCM curves. (Bottom) Time™ 
.series for 1960-2100, relative to 1960, from several 2D-CO2 only 
simulations which u.se time dependent perturbations for the differ- 
ent processes indicated as discussed in Figure Al. The sum of the 
curves in the bottom panel equals the 2D-CO2 only simulation (red 
curve) in the top panel. See text for details. 


Therefore, it appears that the same basic characteristics of the 795 
CCM-simulated BDC and age of air changes are also present 
in the 2D model. 

The CO 2 perturbation shown in Figure 12 (top) mod- 
ifies the model temperature field via the direct IR effect, 
and through changes in surface temperature, latent heating, ooo 


and tropospheric H 2 O which are parameterized from the 
GEOSCCM (Figure Al). These processes are likely due, at 
least in part, to the response of the GEOSCCM hydrological 
cycle to long term SST changes. To examine the relative con- 
tributions of these processes to the 2D model age evolution, 
we ran four additional time dependent COj-only simulations 
in which the surface temperature, latent heating, tropospheric 
H 2 O, and direct IR effect were varied individually, with the 
other CO 2 effects held fixed at 1960 levels. The resulting age 
of air changes are shown in the bottom panel of Figure 1 2, 
relative to 1960 values (we restrict our analysis here to the 
1960-2100 time period when the age changes are largest). 
The increase in latent heating has by far the largest impact, 
contributing ~50% (—0.36 years) of the total C 02 -induced 
mean age change (—0.73 years) from 1960-2100. The im- 
pacts due to tropospheric H 2 O (purple line), surface temper- 
ature (green line), and the direct IR radiative effect (black 
line) are all a factor of '^3 smaller, with age changes of — 0.1 
to —0,14 years from 1960-2100. The magnitude of the di- 
rect IR effect is similar to the GEOSCCM results reported in 
Oman et al. [20091. Note that the sum of the four curves 
in the bottom panel equals the full C 02 -only simulation (red 
curve) in the top panel. 

The ODS, N 2 O, and CH 4 perturbations shown in Figure 
1 2 (top) modify the 2D model temperature field via the small 
direct IR radiative effect (lower atmospheric warming) and 
the larger indirect radiative effects caused by the induced 
ozone and H 2 O changes. These temperature changes modify 
the model zonal wind, wave driving, and BDC, which in turn 
drive the age of air changes shown in Figure 12 (top). 

The impact in the ODS-only simulation (blue solid curve) 
is driven mainly by the strong cooling trend during 1970- 
2000 associated with the ozone hole (Figure B2), which in- 
duces a positive trend in zonal wind (increasing westerlies) 
at high SH latitudes. The accompanying trends (enhance- 
ments) in the planetary wave driving and BDC acceleration 
decrease the age of air. For the ODS-only simulation without 
the ozone hole (blue dashed-dotted curve), the ozone loss and 
cooling is confined to the upper stratosphere globally, with 
only small positive trends in zonal wind and wave driving 
at midlatitudes of both hemispheres, resulting in small BDC 
acceleration and age of air changes in Figure 12. 


6 Photochemical Lifetimes 

Given the significance of the long term stratospheric changes 
caused by GHG and ODS loading, it is useful to examine 
how these changes impact the modeled photochemical life- 
times of certain compounds. The lifetime is important in de- 
termining the length of time over which a molecule of a sub- 
stance will have a significant impact on ozone depletion or 
global warming, and in deriving surface mixing ratio bound- 
ary conditions from emissions estimates for use in atmo- 
spheric models [Kaye et al., 1994]. The lifetime is computed 
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Fig. 13. Time series for 1960-2100 of the photochemical lifetimes 
forNjO.CFC-ll (CFCl3),CFC-12(CF2Cl2),and CCI 4 from the 
base 2D model simulation (black lines). The red lines are from the 
same base simulation except with the loss rates of N 2 O, CFC-1 1 , 
CFC-12, and CCI 4 fixed at seasonally repeating 1960 values. The 
green lines are from a simulation using all chemistry and transport 
fixed at .seasonally repeating 1960 values illustrating the effect of 
the changing atmospheric burden. 


as the atmospheric burden (total number of molecules) di- 
vided by the loss rate, both of which are vertically integrated 820 
and globally/annually averaged. We have recently shown 
the impact of new photolysis cross sections on the model- 
80 S computed lifetime of CCI 4 [Rontu Carton et al., 2010]. Here 
we examine the time dependence of the lifetimes of various 
compounds in more detail. aas 

Figure 13 (black lines) shows the modeled lifetimes of 
N 2 O, CFC-1 1 (CFCI. 3 ), CFC-12 (CF 2 CI 2 ), and CCI 4 from 
810 the base 2D simulation. We restrict our analysis here to the 
1960-2100 time period since the CFC lifetimes are not well 
defined prior to 1960 given that emissions began in the lateeao 
1930s-1940s. While this is not a problem for N 2 O given 
it’s significant natural source, the computed N 2 O lifetime de- 
81 S crease from 1860-1960 is small (143-141 years), consistent 
with the small age of air decrease shown in Figure 12. 


Fig. 14. (Top) The difference between the base and fixed loss rate 
simulations (black lines minus red lines in Figure 13), illustrating 
the lifetime changes due to the changing loss rates. (Bottom) The 
difference between the fixed loss rate and fixed chemistry and trans- 
port simulations (red lines minus green lines in Figure 13), showing 
the lifetime changes due to changes in transport. 


The present day lifetimes shown in Figure 13 (132, 61, 
108, and 5 1 years, respectively, for N 2 O, CFC-1 1 , CFC-12, 
and CCI 4 ), are older than those cited in WMO [2007, 201 1] 
and IPCC [2007]: 1 14, 45, 100, and 26 years. For CCI 4 , the 
older lifetimes shown here do not include soil or ocean toss 
processes. Updated lifetimes for CFC-1 1 and CFC-12 com- 
puted from various models (including the 2D model) were 
presented in Douglass et al. [2008]. That study illustrated a 
strong dependence of lifetime on the modeled circulation and 
age of air, and showed that models with realistic age of air 
simulated a relationship between mean age and the fractional 
release of CFC- 1 1 and CFC- 1 2 that compared well with ob- 
servations. We note that the present 2D model compares well 
with the age of air derived from observations in Figure B4. 

Figure 1 3 (green lines) also shows a simulation in which 
the time dependent loading of N 2 O, CFC-1 1 , CFC-12, and 
CCI 4 impacts only the atmospheric burden used to compute 
the lifetimes, while the model transport and chemistry remain 
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fixed at seasonally repeating 1960 values. This illustrates theaso 
effect of the changing burden on the computed lifetimes. As 
the surface mixing ratio of a substance increases with time, 
an increasing fraction of the total atmospheric burden resides 
in the stratosphere where it is lost, as opposed to the tropo- 
sphere where the loss is zero, and this results in a decrease 
in the computed lifetime, i.e., the lifetime has not reached 
equilibrium with the increasing surface boundary conditions 
fe.g., see Kaye et al., 1994], For the CFCs, emissions be- 
gan in the late 1930s- 1940s and ramped up quickly during®® 
the following decades [Butler et al., 19991, so that there is a 
large impact of the increasing burden on the lifetimes prior to 
~1990. Atmospheric loading of CCI 4 began ~1900 so that 
the influence of the increasing burden on the lifetime is small 
by 1960. This effect on the lifetimes of the CFCs and CCI4 
then levels off as the burdens slowly decrease after 2000 , i.e., 
the lifetime has reached equilibrium with the boundary con- 
ditions. N 2 O has a significant natural source with a slow 
increase throughout 1960-2100 due to anthropogenic activ- 
ity (Figure 1), so that the changing burden has little or no“® 
impact on the computed lifetime throughout the time period. 

The lifetimes are also controlled by the loss rates and the 
rate of transport of a species through the stratospheric loss 
region. The loss rates (photolysis and reaction with O('^D)) 
are impacted by the overhead burden of ozone, and the trans- 
port rates are modified via changes in the BDC as discussed 
previously. Both of these processes incur long term changes 
via the ODS and GFIG loading. 

The red curves in Figure 13 are the same as the base sim- 
ulation except with the loss rates of N 2 O, CFC-1 1 , CFC-12,*”® 
and CCI 4 fixed at seasonally repeating 1960 values. The ef- 
fect of the changing loss rates is then isolated by taking the 
difference between the base and fixed loss rate simulations 
(black lines minus red lines), i.e,, the residual is that due to 
only the changing loss rates. The results are shown relative to®’“ 
1960 in Figure 14 (top) and reveal small impacts for all four 
species. These differences approximately follow the chang- 
ing overhead burden of ozone, with lifetime decreases of ~2 
years (i.e., loss rate increase) for all species from 1960-2000 
as ozone decreases due to ODS loading. As future ozone in-s® 
creases with reduced ODS and increased CO 2 , the impact 
of the loss rates during 2000-2100 increases the lifetimes 
by 6 years for N 2 O, 4 years for CFC-12, and 2 years each 
for CFC-1 1 and CCI 4 . These changes are all <4% over the 
1960-2100 time period. This small impact of the changing^o 
loss rates was also reported in Douglass et al. [2008], 

The effect of the changing BDC is isolated by taking the 
difference between the fixed toss rate and fixed chemistry 
and transport simulations (red lines minus green lines), i.e., 
the residual is that due to only the changing transport. These™ 
results are shown in the bottom panel of Figure 14, again rel- 
ative to I960. The lifetime decreases are largest during 1960- 
2000 when the BDC acceleration is fastest (Figure 12), with 
less pronounced decreases after 2000 as the BDC acceler- 
ates more slowly. The net lifetime reductions for 1960-21 00 9 « 


due to the changing BDC are all 1 1-13%: 16 years for N 2 O, 
1 1 years for CFC-1 1,21 years for CFC-12, and 7 years for 
CCI4. 


7 Summary and Conclusions 

In this paper we use an updated version of our GSFC 2D cou- 
pled model to study long term stratospheric changes caused 
by source gas loading for the 250-year time period, 1 850- 
2100. We use numerous sensitivity simulations in which the 
ODSs, CO 2 , CH 4 , and N 2 O are varied individually to sep- 
arate the relative roles of these gases in driving long term 
changes in ozone, temperature, and age of air. We also 
compare the 2D model with global ozone and temperature 
changes observed during the recent past, and with simu- 
lations from the GEOSCCM to illustrate that the updated 
2D model captures the basic processes that drive long term 
stratospheric changes. A detailed description and evaluation 
of the model are provided in the Appendicies. 

For long term changes in ozone, the impacts due to ODS, 
CO 2 , CH 4 , and N 2 O loading all play important roles in dif- 
ferent regions of the atmosphere. GHG loading becomes 
more important in determining future ozone concentrations 
as the ODS loading diminishes through the 21®‘ century. 
CO 2 cooling increases upper stratospheric ozone via reduc- 
tion in the temperature dependent loss rates. CO 2 loading 
also leads to acceleration of the BDC which redistributes 
ozone in the lower stratosphere, with less (more) ozone in 
the tropics (extratropics). The net CO 2 impact on total ozone 
is a decrease in the tropics and an increase in the global av- 
erage and at midlatitudes, with a larger enhancement in the 
NH. For 1850-2100, CO 2 loading has the largest individual 
impact on global ozone in the upper stratosphere and on the 
total column, causing increases of 20% and 4.2%, respec- 
tively. N 2 O loading and the subsequent increase in the NOx 
ozone loss maximizes near 35 km (^^0 .5%/decade), with a 
global total column decrease of —2.7% from 1850-2100. 

Methane loading impacts ozone via several processes, and 
we illustrate the sensitivity of these processes in determin- 
ing the net result on ozone. The simulations, which use fixed 
tropospheric NOy, reveal that CH 4 loading leads to a net in- 
crease in total ozone at all latitudes throughout 1850-2100. 
Global total ozone increases by 2-2.5% from 1850 to the 
latter half of the 21®‘ century due to CH 4 loading. How- 
ever when coupled with time dependent ODS changes, the 
methane impact on ozone for present day conditions is nearly 
a factor of two larger (-F 6 DU) compared to the 4-3 DU re- 
.sponse calculated with very low ( 1 850) levels of Cly . 

In the lower mesosphere, the CH 4 -induced HOx-ozone 
loss dominates the ozone time dependence. Our simulations 
reveal that this process resulted in significant ozone reduc- 
tions prior to 1960, so that one-third of the total decrease in 
ozone during 1850-2050 occurred by 1960. 
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In the stratosphere prior to 1960, the simulated ozone 
changes are relatively small, but are not insignificant. ODS 
emissions, primarily due to carbon tetrachloride, cause 1% 
depletion in global total ozone from 1 860- 1 960, with the vast 
945 majority of this occuring after ~1920. From 1860-1960, in-iooo 
creasing CO 2 loading results in ozone increases of 2 % in 
the upper stratosphere and 0.5% in the total column, with 
CH4 and N 2 O causing total ozone changes of +0.75% and 
—0.5% respectively, during this period. The net result is a 
950 broad ozone maximum during the 1 920s- 1 930s, as the CO2-1005 
and CH4- induced ozone increases outpace the ozone losses 
caused by N 2 O and CCI4 emissions. This proceeds the de- 
cline in ozone driven mainly by ODS loading, which be- 
comes much more rapid after ~ 1970. 

955 We also present a quantitative analysis of the steady stateimo 
impacts of GHGs on ozone for 1850, 1950, 2000, and 2100 
atmospheric conditions. These calculations have a significant 
dependence on the background conditions, reflecting satura- 
tion effects as well as the chemical interactions of the pertur- 
960 bations. Using the ozone depletion potential (ODP) concept 
for GHGs, we compute present day ODP values of +0.019 
(N 2 O), -0.00065 (CH 4 ), and -0.0039 (CO 2 ), relative to 
CFC-1 1. To more easily interpret the relative impacts of theois 
GHGs on ozone, we also derive steady state sensitivity fac- 
965 tors in terms of the mixing ratio change of each GHG. These 
calculations reveal that based on the AIB scenario, CH 4 had 
the largest impact on ozone for year 2000 conditions, owing 
to it’s strong coupling with Cly . However, CO 2 is projectedio 2 o 
to have the largest impact by 2100 , and will have twice the 
970 magnitude of the N 2 O impact. 

The simulated changes in stratospheric temperature for 
1850-2100 are mostly controlled by CO 2 cooling, with the 
reduced ozone heating caused by ODS loading also playingio^s 
an important role from ~1980 through the first half of the 
975 2T’* century. The impact of CH 4 and N 2 O loading are rel- 

atively small below ~45 km. However, the cooling due to 
CH 4 -induced H 2 O increases becomes significant above ~50 
km, and the resultant temperature changes are comparable toioso 
those induced by the ODS and CO 2 loading. 

980 In the 2D model, the long term changes in surface tem- 
perature, latent heating, and tropospheric H 2 O are param- 
eterized based on the GEOSCCM simulations and the CO 2 
loading. These processes are likely due, at least in part, to theoss 
response of the GEOSCCM hydrological cycle to long term 
985 SST changes. Parameterization of these processes enables 
the 2D model to simulate long term changes in the BDC and 
age of air which are consistent with the GEOSCCM. Changes 
in GHG and ODS loading ultimately impact the BDC and ageio 4 o 
of air through changes in the radiative and temperature dis- 
990 tributions. We estimate that changes in ODS concentrations 
account for 65% of the mean age decrease from 1980-2005, 
with 25% due to CO 2 loading, and 10% due to CH 4 and 
N2O. Our simulations also indicate that the BDC and ageuMs 
changes caused by ODS loading are mainly due to the forma- 
995 tion of the ozone hole. CO2 loading becomes the dominant 


source of mean age change after ■'^2025, and explains ~ 75 % 
of the total age reduction ( 1 year) from 1 860-2 1 00 , with 12 % 
and 5% explained by N 2 O and CH 4 loading, respectively. 

We also examined the time dependent photochemical life- 
times of N 2 O, CFCI 3 , CF 2 CI 2 , and CCI 4 , and found that 
the impact of the BDC acceleration is significant, and causes 
the lifetimes to decrease by 11-13% from 1960-2100. The 
impact of the changing loss rates is generally small and fol- 
lows the time dependent changes in stratospheric ozone via 
the changes in photolysis and 0(^D). This effect decreases 
the lifetimes by ^1 years from 1960-2000, and increases the 
lifetimes by 2-6 years (3-4%) from 2000-2100, The model 
calculations also allow us to separate these geophysical im- 
pacts on the lifetimes from the artifacts caused by the dise- 
quilibrium between the rapidly increasing atmospheric bur- 
den and the stratospheric loss. 


Appendix A 

GSFC Coupled 2D Model 

The GSFC 2D coupled chemistry-radiation-dynamics model 
was originally discussed in Bacmeister et al., [19951 and has 
been frequently used in stratospheric assessments [WMO, 
20071, and studies pertaining to the chemistry-climate cou- 
pling of the middle atmosphere fe.g., Rosenfield et al., 1997; 
Rosenfleld et al., 20021. While this model was not included 
in the recent CCMVal activity [SPARC CCMVal, 20101, sev- 
eral of the model components are very similar to those used 
in the GEOSCCM which was evaluated in CCMVal. These 
components include: the infrared radiative transfer scheme 
[Chou et al., 20011; the photolytic calculations [Anderson 
and Lloyd, 1990; Jackman et al., 19961; and the microphys- 
ical model for PSC formation [Considine et al., 19941. As 
discussed in Fleming et al. [20071, the model now uses 
an upgraded chemistry solver that computes a full diurnal 
cycle for 35 fast chemical constituents. This scheme was 
shown to be in good agreement with photochemical steady 
state box model calculations [Park et al., 19991. The latest 
Jet Propulsion Laboratory (JPL) recommendations are used 
for the photolytic cross sections and reaction rate constants 
[Sander et al„ 2006[. 

The model domain extends from the ground to approxi- 
mately 92 km. The chemistry calculations are done on a grid 
resolution of 4" latitude by 1 km altitude. We have found that 
for most applications, the model radiation and dynamics cal- 
culations can be adequately done on a somewhat coarser grid 
of ~ 4.9" latitude by 2 km altitude. Using a finer resolution 
only adds to the computational burden but does not improve 
the model dynamical simulations. 

We have recently made extensive upgrades to the model 
solar radiation and dynamical modules, which are described 
in the following sections. An evaluation of the model tem- 
perature and transport are then provided in Appendix B. 
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Al Wave Parameterizations two 

The horizontal mixing (Kyy) and momentum deposition due 
1050 to dissipating planetary waves is computed using a linearized 
parameterization similar to that described by Garcia [19911. 
Previously, the parameterization used a surface boundarynos 
condition based on a representation of the topographic forc- 
ing of planetary waves. It was necessary to include adjustable 
1055 amplitude efficiency factors for each wave number to obtain 
reasonable seasonal variations of the zonal winds and chemi- 
cal fields [Rosenfield et al., 19971, In the present model, this 
lower boundary condition is based on a geopotential height 
climatology for 750 mbar as a function of latitude and sea-^^^ 
1060 son, derived from the National Centers for Environmental 
Prediction (NCEP) reanalysis-2 data [Kistler et al., 20011 
averaged over 1979-2007. This provides a more complete 
representation of the lower boundary forcing of planetary 
waves, e.g., land-ocean contrasts, in addition to the topo- 
1065 graphic forcing. We solve for planetary zonal wave numbers 
1-4. Compared to that obtained previously, the new method- 
ology provides more realistic model simulations of planetary 
wave drag and mixing with no artificial wave amplitude ad- 
justment factors necessary. 

'' 1120 
1070 The model also includes the off-diagonal eddy diffusion 

of constituents (Ky^), following the methodology used in our 
GSFC 2D fixed transport model. This follows the assumption 
that horizontal eddy mixing is directed along the zonal mean 
isentropes, and projects the Kyy mixing rates onto isentropic^,^^ 
1075 surfaces [Plumb and Mahlman, 1987; Newman et al., 19881. 

The momentum deposition due to gravity wave breaking in 
the mesosphere is parameterized using a ray tracing calcula- 
tion for waves with non-zero phase speeds, and a cubic drag 
law for zero-phase speed mountain waves [Bacmeister et al.,,^,, 
1080 19951. This enables the gravity wave momentum flux to be 

interactive with the evolution of the zonal mean flow. How- 
ever to obtain proper tracer simulations, we found it neces- 
sary to specify the model vertical eddy diffusion rates 
which are taken from the GSFC 2D fixed transport model,, 3 ^ 
1085 as a function of latitude, height, and season [Fleming et al., 
20071, For the upper stratosphere and mesosphere, these are 
based on the gravity wave parameterization originally devel- 
oped by Lindzen [19811 and modified by Holton and Zhu 
[19841. For the troposphere and lower stratosphere, is 
1090 based on the zonal mean temperature lapse rate as computed 
from a multi-year average of NCEP reanalysis-2 data. hm 

A2 Radiative TVansfer 

For the absorption of solar radiation in the ultraviolet and vis- 
ible, we now compute the heating rates consistently with thaws 
1095 model incident solar flux and photolysis calculations [Jack- 
man et al., 19961. These are computed over a full diur- 
nal cycle with a much finer spectral resolution compared to 
the broad band calculations used previously [Strobel, 1978; 
Chou and Suarez, 19991. These heating rate calculationsnso 


have the further advantage of utilizing the most current rec- 
ommendations for the photoly tic cross sections [Sander et al ., 
20061 . This new methodology also gives model temperatures 
that are in somewhat better agreement with observations. 

For ozone, the heating rate in this way is computed by [e.g, 
Brasseurand Solomon, 19861: 

dT 1 f 

— = — [03]y^(T(O3)F,,AdA (Al) 

where p is the total atmospheric density, Cp is the specific 
heat of dry air at constant pressure, [O 3 ] is the ozone number 
density, A is the wavelength, cT(_Oz) *s the ozone absorption 
cross section, and Fg,x is the incident solar flux (enhanced 
or reduced) as a function of wavelength at each model grid 
point. In addition to ozone heating, we also include the heat- 
ing due to absorption by O 2 , which is important in the meso- 
sphere, and NO 2 which is of secondary importance in the 
middle stratosphere. For these calculations, we assume that 
all of the solar radiation absorbed is immediately realized 
as thermal energy. This is a good approximation below ~80 
km, i.e., the region of interest for the current study, where the 
chemical recombination of 0(®P) is very fast [e.g, Brasseur 
and Solomon, 19861. 

We also include the minor absorption of solar radiation 
by ozone in the infrared, and by water vapor in the in- 
frared and visible based on the parameterization of Chou and 
Suarez [19991. For these calculations, we use surface reflec- 
tivity values as a function of latitude and season based on 
the TOMS climatology compiled by Herman and Celarier, 
[19971. 

For the thermal infrared radiative transfer, we use the pa- 
rameterization of Chou et al. [2001], which is the same as 
that used in the GEOSCCM. This includes the contributions 
due to O 3 , CO 2 , H 2 O, CH 4 , N 2 O, CFC-11, CFC-12, and 
HCFC-22. For both the solar and thermal IR calculations, 
the 2D model includes zonally averaged cloud parameters 
based on a multi-year average of output from the Whole At- 
mosphere Community Climate Model version 3 (WACCM3) 
[e.g., Garcia et al., 20071. 

A3 Model Treatment of Longitudinal Variations 

Previously, the model-generated zonal mean temperatures 
were used to compute the gas phase and heterogeneous reac- 
tion rates. In the new model version, reaction rates are now 
computed using a longitudinal temperature probability dis- 
tribution which is generated from the model-computed plan- 
etary wave fields for zonal waves 1-4 (section Al). The rates 
for the reactions at each model grid point are computed once 
per day by summing the rates computed for each temperature 
in the distribution weighted by the probability of occurrence 
of that temperature. Using the temperature probability dis- 
tribution instead of the zonal mean temperature is especially 
important for the heterogeneous chemical reactions as these 
can have significant non-linearities in temperature. 
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For the calculation of polar stratospheric cloud (PSC) for- 
mation, we utilize the parameterization described in Consi- 
dine et al. [19941, using longitudinal temperature probability 
distributions derived from the NCEP reanalysis-2 data aver- 
aged over 1979-2006. This climatological average distribu- 
tion is used for all years in the simulations. This method- 
ology does not allow for the interaction between PSC for- 
mation and the chemical/dynamical time evolution of the 
model stratosphere. However given the strong temperature 
non-linearity of PSC formation, we found it necessary to use 
the observed temperature distributions rather than the model 
temperatures to properly simulate PSC formation. 

A4 Tropospheric Parameterizations 

Accounting for the hydrological cycle and surface boundary 
layer processes is important to properly simulate the dynam- 
ical and chemical distributions of the troposphere and lower 
stratosphere. Since the 2D model framework is inadequate to 
simulate most tropospheric processes interactively, we spec- 
ify the surface temperature, tropospheric water vapor and la- 
tent heating. As described in the following, we first generate 
monthly and zonal mean climatologies of these parameters, 
and then add on long term changes parameterized in terms of 
the atmospheric CO 2 loading. 

The surface temperature seasonal cycle as a function of lat- 
itude is based on the NCEP reanalysis-2 data averaged over 
1979-2006. Tropospheric latent heating as a function of lati- 
tude, height, and season is based on a multi-year average of 
output from WACCM3 simulations. The model water vapor 
seasonal cycle in the upper troposphere (12-16 km) is based 
on the UARS reference atmosphere (UARSRA) compiled by 
Randel et al., [20011. Below 12 km, water vapor is derived 
from a 2 1 year average (1981 -200 1 ) of relative humidity data 
from the European Center for Medium Range Weather Fore- 
casts (ECMWF) updated reanalyses (ERA-40). Values from 
the UARSRA and ERA-40 data sets are functions of latitude, 
height, and season and are blended over several pressure lev- 
els to obtain a .smooth transition in the vertical. Water va- 
por everywhere above the tropopause is computed in the 2D 
model (see Figure B6). 

In addition to the seasonal variations, the surface tem- 
perature, tropospheric water vapor and latent heating un- 
dergo substantial long term changes as simulated by the 
GEOSCCM. This is illustrated in Figure Al which showsizos 
zonally averaged time series of GEOSCCM simulations for 
1950-2100 (black curves) at the locations indicated. These 
time series have been deseasonalized and smoothed to re- 
duce the interannual variability of the GEOSCCM. These 
long term changes are highly correlated with the time de-i 2 io 
pendent surface boundary condition of CO 2 (bottom panel) 
and are likely a response to the warming of the troposphere 
and sea surface temperatures caused by the increased atmo- 
spheric CO 2 loading. To represent these long term changes 
in the 2D model, we compute a sensitivity factor to the CO 21215 
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Fig. Al. Time series for 1950-2100 of zonally averaged and desea- 
sonalized surface temperature, water vapor mixing ratio, and latent 
heating from the GEOSCCM (black curves), and the corresponding 
fits to the surface CO 2 boundary condition (red curves) for the lo- 
cations indicated. Inclusion of the seasonal cycle to the CO 2 fits is 
depicted by the orange shading. The surface CO 2 boundary condi- 
tion is shown in the bottom panel. See text for details. 


boundary condition for each parameter at each model grid 
point. The time dependent value of each parameter is then 
determined by the CO 2 concentration at each time step mul- 
tipled by the sensitivity factor. As depicted by the red curves 
in Figure Al , this smoothly-varying fit to the CO 2 boundary 
condition allows us to remove the unwanted artifacts of the 
GEOSCCM interannual variability. This also allows us to 
extrapolate these quantities to years prior to 1 950 based on 
the CO 2 loading, assuming the same sensitivity to CO 2 as 
for 1950-2100. The orange shading in Figure Al shows the 
climatological seasonal cycle for each parameter added onto 
the long term variation for input into the 2D model. 
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Fig. Bl. Latitude-height cross-sections of the February monthly 
mean temperature (K) from the MERRA data (top), 2D model sim- 
ulation (middle), and the difference , model minus MERRRA (bot- 
tom panel), all averaged over 1979-2009. The contour intervals are 
lOK for the top and middle panels and ±2 K for the bottom panel. 


Appendix B 

2D Model Temperature and Tracer Comparisons 

In this Appendix, we provide an evaluation of the 2D model 
1220 transport fields by comparing the simulations with observa- 
tions of several tracers including ozone. We also evaluate the 
model temperature simulations (climatology and trends). 1235 

Bl Temperature 

The model temperature field for February is shown in Figure 
1225 Bl, along with the MERRA meteorological analyses (sec -1240 
tion 4), both averaged over 1979-2009. The largest model 
differences (bottom panel) occur at high SH latitudes in the 
upper troposphere (— lOK), and in the upper stratosphere at 
high NH latitudes (-1-I4K). However overall, the model is in 
1230 reasonable agreement with the MERRA data, as the modeli245 
differences are mostly within ±5K. 


Fig. B2. (Top panels) Annually averaged temperature trend for 
2000- 1 960 from the GEOSCCM (left) and 2D model (right) base 
simulations (all source gases varied time dependently). The trends 
are derived from linear least squares fits to the annually averaged 
time series. (Bottom panels) As in the top panels except for the 
2095-2005 temperature difference (K/decade) using 10-year aver- 
ages centered on 2095 and 2005 to reduce the effects of interan- 
nual dynamical variability in the GEOSCCM. The contour interval 
is ±0.1 K/decade. 


The past and future temperature changes from the 2D 
model and GEOSCCM are shown in Figure B2. The 2D 
model captures most of the latitude-height variations sim- 
ulated in the GEOSCCM. The main discrepancies occur at 
high latitudes where the 2D model somewhat underestimates 
the large temperature changes simulated by the GEOSCCM 
associated with the ozone hole, i.e., past cooling and fu- 
ture warming. Also, the GEOSCCM simulates a mid-upper 
stratospheric warming at high SH latitudes for 1960-2000, 
which was shown to be a dynamical response to the ozone 
hole fStolarski et al., 2006). The 2D model simulates this 
feature only very weakly (top right). These discrepancies are 
likely due to the 2D model not fully resolving the large zonal 
asymmetries characteristic of the polar region, as well as the 
known high ozone bias at high latitudes in the GEOSCCM 
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Fig. B3. As in Figure B2, except for the global (90°S-90“N) av- 
erage profiles. The top panel includes the trend derived from ra- 
diosonde data (RATPAC-A) [Free et al., 2005]. 


[Pawson et al., 2008]. 

The corresponding global average vertical profiles (Figure 
B3) also show good agreement between the 2D model and 
1250 GEOSCCM as well as radiosonde data for 1960-2000 from 
the Radiosonde Atmospheric Temperature Products for As- 
sessing Climate (RATPAC-A) [Free et al., 2005]. We note 
that the global average GEOSCCM stratospheric tempera- 
ture trends were also found to be in reasonably good agree- 
1256 ment with those derived from SSU and MSU satellite data 
for 1979-1999 [Stolarski et al., 2010]. 

In the troposphere, the 2D model simulates warming 
throughout 1960-2100, which is due mainly to the param- 
eterized long term changes in surface temperature and latenhsrs 
Km heating shown in Figure Al . The magnitude of the warming 
in the tropics is somewhat underestimated compared with the 
GEOSCCM in both the past and future (Figure B2), although 
the 2D-simulated warming over 1 960-2000 compares favor- 
ably with the radiosonde data in the global average in Figuroaso 
1265 B3. 


Fig. B4. Age of air at 20 km derived from ER-2 aircraft measure- 
ments of SFe (asterisks) and CO 2 (triangles), and vertical profiles 
of the age of air derived from balloon measurements of SFe (aster- 
isks, plus signs) and CO 2 (triangles) at the latitudes indicated. Ages 
derived from these measurements have been adapted from Hall et al. 
(1999). Also shown are simulations from the 2D model (red line). 
The age is taken relative to the tropical tropopause. 


B2 Age of Air 

Stratospheric age of air is a widely used diagnostic that tests 
the overall fidelity of model transport [e.g.. Hall et al., 1999]. 
Figure B4 shows the mean age of air at 20 km derived from 
aircraft measurements of SFe (asterisks) and CO 2 (trian- 
gles), and a series of vertical profile measurements of SFg 
and CO 2 made from balloon flights in three latitudes zones 
[e.g., see Hall et al., 1999], We note that differences in the 
observations at the middle and higher latitudes may reflect 
photo-chemical influences on SFg which would cause an 
overestimation in the inferred ages [Hall and Waugh, 1998]. 
Some of the older age measurements at 65“N may also reflect 
remnants of the winter polar vortex [Ray et al., 1999]. 

Figure B4 also shows the age of air derived from the 2D 
model simulation averaged over the 1990s (red line). The 
age of air in the model is computed from a ’’clock” tracer 


Altitude (km) Altitude (km) 
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Fig. B5. Latitude-height cross-sections of September N 2 O aver- 
aged over 2004-2009 from AURA/MLS (top) and the 2D model 
(bottom). Contour interval is 20 ppbv and includes the 10 ppbv 
contour. 


that has a surface boundary condition linearly increasing with 
time, with no other chemical production or loss. This is es- 
sentially identical to the age obtained from simulations of 

1285 SFe or CO 2 as done, for example, in Hall et al. [19991.'«i 
The model somewhat underestimates the observed age in the 
mid-latitude Northern Hemisphere (NH) (see the 20 km and 
40°N panels), and at the high latitudes of the NH above 30 
km. However for the most part, the model simulates the ab- 

12 S 0 solute values and the latitudinal and vertical gradients of theisos 
observations fairly well. This illustrates that the model trans- 
port rates in the stratosphere, i.e., the relative magnitudes of 
vertical motion and horizontal mixing, are generally realistic. 

B3 NaOandHaO 

1295 As further evaluation of the model transport, we show 
latitude-height cross sections of September N 2 O from the 
model and the AURA/MLS data averaged over 2004-2009 
(Figure B5), and March H 2 O from the model and thesis 
UARS/HALOE data both averaged over 1994-2004 (Figure 


Fig. B6. Latitude-height cross-sections of March H 2 O averaged 
over 1 994-2004 from UARS/HALOE (top) and the 2D model (bot- 
tom). The top panel includes data from AURA/MLS in the polar 
regions (averaged over 2004-2009) where HALOE lacks data cov- 
erage. Contour interval is 0.2 ppmv. The red dashed line separates 
the regions where the model H 2 O is computed (above) and pre- 
scribed to the HALOE climatology (below). 


B6). Figure B6 includes data from AURA/MLS in the polar 
regions (averaged over 2004-2009) where HALOE lacks data 
coverage. The model shows good overall agreement with the 
data in reproducing transport sensitive features in the merid- 
ional plane, including the horizontal and vertical gradients. 
For example, the model qualitatively simulates the region 
of strong horizontal mixing during late winter/early spring 
at midlatitudes of both hemispheres. This is especially pro- 
nounced in the SH during September at 20-40 km in the N 2 O 
field. The model tends to underestimate this mixing in the 
SH mid-upper stratosphere, as is also seen in the midlatitude 
vertical profile in Figure B7 (top). Here the model compares 
well with the MLS N 2 O below ~27 km, but underestimates 
the data above this level, which is due to weaker than ob- 
served mixing of high-N 20 air from lower latitudes. 

In the Antarctic lower stratosphere, the observed N 2 O pro- 
file in September (Figure B7 , bottom) is mainly a result of de- 
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H2O+2CHJ seasonal cycle Equator 
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Fig. B7. Vertical profiles of September N 2 O at 30°S-50"S (top) 
and 80°S-88°S (bottom) from AURA/MLS observations (black tri- 
angles) and the 2D model (red curve) averaged over 2004-2009. 
The error bars (Icr) denote the combined effects of measurement 
uncertainty and interannual variability. The Antarctic MLS profile 
is adapted from SPARC CCMVal 1 20 101. 


Fig. B8. Equatorial profiles of the amplitude and phase lag of the 
.seasonal cyle in the quantity H 2 O-I- 2 CH 4 from UARS/HALOE 
(black asterisks) and the 2D model (red). Amplitudes are relative 
to the values at 16.5 km, and the phase lag is defined to be zero at 
16.5 km. Values are averaged over 1994-2004 for both the data and 
model. 


scent within the vortex occurring throughout the winter, withsss 
little influence of air in-mixing from midlatitudes fSPARC 
CCMVal, 2010]. The good model-data agreement here il- 
1320 lustrates that the magnitudes of vortex descent and isolation 
from midlatitudes in the simulation are generally realistic. 

The HALOE H 2 O data in Figure B 6 indicate strong pole-i 3 « 
ward and downward transport of very dry air from the tropics 
to midlatitudes just above the tropopause fe.g., Randel et al., 
1325 2001]. The model H 2 O is set to the HALOE climatology 
in the upper troposphere, below the red dashed line in Fig- 
ure B 6 . Above this level, the H 2 O field is computed in tha 34 s 
model, and reveals that the model transport is resolving fairly 
well this strong poleward and downward transport from the 
1330 tropical tropopause region. 

In Figures B5 and B 6 , the model also resolves the isola- 
tion of the tropics in the lower stratosphere, as indicated byiaso 
strong horizontal gradients in the subtropics, and a region of 
low water vapor concentrations at 20-30 km over the equator 


associated with the “tape recorder” signal [e.g., Mote et al., 
19961. This feature reflects the slow upward propagation of 
the water vapor seasonal cycle from the tropical tropopause, 
and simulation of this feature provides a good diagnostic of 
model transport. 

Figure B 8 shows the amplitude variation and phase lag 
versus altitude of the seasonal cycle in H20-f 2CH4 at the 
equator relative to the tropopause from HALOE data (black 
asterisks). This quantity is quasi-conserved and accounts 
for both the H 2 O seasonal cycle propagation and the slow 
photochemical conversion of CH 4 into H 2 O in the strato- 
sphere. The amplitude attenuation and phase lag with in- 
creasing height reflect the strength in the upwelling of the 
Brewer-Dobson circulation (BDC) combined with the rate of 
vertical diffusion and entrainment of air from mid-latitudes 
[Hall et al., 1999L The model (red line) shows an increas- 
ingly longer phase lag compared with the data above ~27 
km, possibly reflecting weaker BDC upwelling in the trop- 
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Fig. B9. Latitude-height cross-sections of annually averaged ozone 
expressed in Dobson Units per kilometer. The observations (top 
panel) are from the climatology compiled by McPeters et al. [20071 
covering the time period 1988-2002. Also shown are the 2D model 
simulation averaged over 1 988-2002 (middle panel), and the differ- 
ence, model minus observations (bottom panel). For the top and 
middle panels, the contour interval is 2 DU/km and includes the 
3 contour level. For the bottom panel, the contour interval is ±1 
DU/km. 


ical middle stratosphere than indicated in the observations. 
However overall, the model shows mostly good agreement 

1355 with the HALOE data in simulating this seasonal cycle prop- 
agation. This, combined with the good agreement in the trop- 
ical age profile (Figure B4), suggests that the model transport^^^j 
rates in the tropical lower-middle stratosphere appear to be 
fairly realistic. 

1360 B4 Ozone 

1380 

Figure B9 shows latitude-height cross sections of annu- 
ally averaged ozone from an observational climatology (top 
panel), the model (middle panel), and the difference, model 
minus data (bottom panel) expressed in DU per kilometer. 

1365 This unit is proportional to the number density per cm^ di-isBs 
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Fig. BIO. Month-latitude cross-sections of total column ozone av- 
eraged over 1988-2002 from ground-based observations (top panel, 
updated from Fioletov et al. [2(X)2]) and 2D model simulation (bot- 
tom panel). The contour interval is 20 DU. 

vided by a constant, and is a direct measure of the contri- 
bution versus altitude to the total column. The climatol- 
ogy is ba.sed on a combination of ground based and satel- 
lite data covering the time period 1988-2002, as compiled by 
McPeters et al. [2007]. The model results are also averaged 
over 1988-2002. 

The simulation compares relatively well with the data in 
most regards. The model qualitatively reproduces the ob- 
served latitudinal and vertical gradients in most places, as 
well as the magnitude of the ozone amounts. There are some 
regions of small discrepancy; for example at 15-20 km in 
the SH midlatitudes where the model slightly overestimates 
ozone, and in the tropics above 30 km associated with the 
ozone deficit region where the model underestimates ozone 
[e.g., Jackman et al., 1996]. The model also underestimates 
ozone at mid-high NH latitudes near 10-15 km which is 
likely due to excessive horizontal mixing in this region, i.e., 
in-mixing of low-ozone air from the tropical troposphere. 
However, these model-measurement differences are now sig- 
nificantly smaller compared with previous model versions. 
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The corresponding season-latitude sections of total col-1440 
umn ozone also averaged over 1988-2002 are shown in Fig- 
ure BIO. The observations are from ground-based measure- 
ments updated from Fioletov et al. [20021. Again the model 
1390 shows good overall agreement with the data in reproducing 
the absolute total ozone magnitudes, the latitudinal gradients,"”* 
and the seasonal variations. The main discrepancy occurs at 
NH high latitudes where the model tends to underestimate 
total ozone throughout the year by ~20-50 DU. This reflects 
1396 the model underestimation of ozone at high NH latitudes in, 45 „ 
the 10-15 km region seen in Figure B9 (bottom panel). 
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